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Stanene, composed of tin atoms arranged in a single layer, is the tin analogue of 
graphene and past studies predicted it to be a topological insulator. An energy band 
gap (of 0.1 eV) was obtained in previous calculations for the buckled honeycomb 
structure of stanene and, thus, phonon-mediated superconductivity in this material 
is ruled out. In this work we investigated, from first principles calculations within 
density functional theory (DFT), the possibility of producing phonon-mediated su¬ 
perconductivity in stanene by doping the material. It was found that doping with 
calcium (lithium) leads to superconductivity, albeit, with a very low superconducting 
transition temperature Tc of ~ 0.7 K (~ 1.3 K), even lower than the value (3.7 K) 
for bulk /3-Tin. 
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I. INTRODUCTION 

Stanene is a two dimensional sheet composed of tin (Sn) atoms arranged in a honeycomb 
lattice, similar to graphene. Ab initio calculations^ indicate that, unlike graphene, a low- 
buckled conhguration of stanene is more stable compared to the planar geometry. This is as 
a result of the weak vr — tt bonding between the atoms. The buckling enhances the overlap 
between vr and a orbitals and stabilizes the system.^ Theoretical studies indicate that free 
standing stanene is a zero band gap semiconductor in the absence of spin orbit coupling 
(SOC). The inclusion of SOC leads to a band gap of about O.leV.^’^ Chemical functional¬ 
ization has been predicted^’^ to be capable of increasing the band gap to up to 0.35 eV. 
With this capability of tuning its bandgap, stanene could be useful as a semiconductor in 
low-power-consumption electronics. 

From the predicted electronic zero/non-zero band gap of stanene, one can conclude that 
phonon-mediated superconductivity is not possible in this material, just like its carbon 
counterpart graphene is not superconducting. Recently, some calculations^ have, however, 
shown that superconductivity can be induced in graphene by n-doping at a carrier density 
exceeding 10^® cm“^. Superconductivity has been produced in graphene by doping its surface 
with alkaline metal ad atoms.®’® This is reminiscent of the induction of superconductivity 
in graphite intercalated compounds (GICs).® Specihcally, it was found that lithium-covered 
graphene has a higher superconducting transition temperature Tc than its calcium-covered 
counterpart unlike the situation in GICs where the opposite is the case.® Yang et al. showed 
that the superconductivity in GIGs is as a result of the interaction between the interlayer (IL) 
states and the tt* states, and that, if an ad atom superlattice could be created with graphene, 
a similar situation could lead to superconductivity.^ This was conhrmed experimentally for 
doped graphene in Ref. 5. 

With the electronic similarities between graphene and stanene, namely, the existence of 
a Dirac cone at the so called K point and the presence of a tunable band gap^’®, could it 
be possible to produce superconductivity in stanene in a manner similar to the induction of 
superconductivity in graphene? Doped stanene might superconduct at an even higher tran¬ 
sition temperature Tc than doped graphene owing to the fact that bulk tin (white Tin) is su¬ 
perconducting at low temperatures without doping, in contrast to bulk graphene (graphite) 
which is not superconducting. Graphite can mainly be made superconducting by creating 
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interlayer states in it via metal intercalation. 

To this end we examine, from first principles calculations, within Density Functional 
Theory (DFT), the possibility of inducing superconductivity in stanene by creating IL states 
in it by doping. We doped the surface of stanene with an alkaline and an alkaline earth metal, 
lithium and calcium, respectively and predicted Tc for these doped systems. We note that 
decorated stanene (SnX, with X= OH, I, Cl, Br) have been predicted to be topological 
superconductors when doped with silver atoms^ but the superconducting property of SnX 
itself was not predicted but only assumed in that study. 

Stanene has recently been synthesized.® It was fabricated on a Bi 2 Te 3 (lll) substrate by 
the process of molecular beam epitaxy.® This experimental realization of stanene should 
open up the possibility of verifying the different predictions made for the material. 


II. METHODS 

We first determined the optimum structure of stanene and computed its electronic band 
structure. The calculations were performed using Density Functional Theory (DFT) with the 
gradient corrected exchange correlation functional of Perdew, Burke, and Ernzerhof (PBE).^® 
The QUANTUM ESPRESSO^^ suite of codes was used. Plane waves with a kinetic energy 
cutoff of 55 Ry were used in expanding the wavefunctions and core electrons were excluded 
through the use of a norm-conserving pseudopotential. We inserted a vacuum layer with 
a width larger than 20A in the non-periodic direction of the system. Since the material 
under consideration is a semimetal, the electronic properties were determined by integration 
in the first Brillouin zone (BZ) using a Gaussian smearing having a width of 0.05 Ry. 
All conhgurations and their respective atomic positions were optimized ensuring that the 
maximum force on any atom was smaller than 10“^eV/A. The electronic band structure 
was calculated along the path connecting the special high symmetry k-points F, M, K in the 
first BZ as shown in Sect. Ill below. The electronic density of states (EDOS) was computed 
using a hner mesh of the first BZ with electron momentum of 27 x 27 x 1 and 21 x 21 x 1 
for stanene, and doped stanene respectively. 

Next, the phonon band structure and phonon density of states (PDOS) were calculated 
in the framework of density functional perturbation theory (DFPT) as implemented in the 
QUANTUM ESPRESSO codes.These calculations showed, among other things, the sta- 
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bility of the buckled stanene structure. The electron-phonon coupling (EPC) constants A 
were also determined but only for doped stanene since undoped stanene does not super¬ 
conduct. The necessary summations over the electron momenta and phonon momenta (see 
Eqs. 4 and 5 below) in computing the EPC constants were carried out using an electron 
momentum mesh of 32 x 32 x 1 and a phonon momentum one of 8 x 8 x 1. 

We used the following formula by Allen and Dynes^'^ to determine the superconducting 
transition temperature T^. 


T 


Ulog f 1.04(1 -l- A) A 

(,-A-;.*(l + 0 . 62 A)J’ 


( 1 ) 


where, uiog, A and /i* are the logarithmic average of the phonon energy, electron-phonon 
coupling constant and the electron-electron Coulomb repulsion parameter, respectively. The 
hrst two quantities were calculated as: 


Ulog = exp 




a^F{u) Incu—^ 


( 2 ) 


^oo 2 


A = 2 


a^{u)F{u) 


dw 


u 


where, the Eliashberg function a^F is given by 


(3) 


“ jvmw iv, A] (4) 

^ ^ k,k',D 

The terms that appear in the expression are dehned as follows: A^(0) is the electronic 
density of states per spin at the Fermi level, i.e., A^(0) = <5(ek) and A^k and A"q represent 

the total number of k and q points respectively. The Kohn-Sham eigenvalue with reference 
to the Fermi level is and is the electron-phonon coupling matrix element. The Dirac 
delta distributions 5(ek) and (5(ek') restrict the electrons to the Fermi surface, Uc^^u is the 
phonon frequency at the phonon wave vector q of the phonon branch of index z/. Here, k 
and k' indicate both the electron wave vector and band index, and spin index, q = k' — k 
indicates both the phonon wave vector and band index. In our calculations the Dirac delta 
distribution of Eq. 4 above was replaced by a Gaussian broadening function having a width 
of 0.01 Ry. The variable is the coupling strength and F{u) is the phonon density of state. 
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i.e., 

q 

III. RESULTS AND DISCUSSION 

The optimized lattice parameters are shown in Table I along with the results obtained 
from calculations in the Literature. The variable a is the lattice constant, d is the out-of- 
plane buckled height, and agn-sn is the Sn-Sn bond length. The values calculated in this 
work are consistent with other theoretical predictions in the literature. We also found the 
buckled stanene structure to be stable. The computed phonon band structure of pristine 
stanene had no imaginary frequency. This stability is in agreement with Ref. 2,3 but differs 
from the results presented in the supplementary information of Ref. 15. 

Next, we present the electronic band structure and the electronic density of states of 
free standing stanene. The band structure and the density of states in Fig. 2 are in good 
agreement with the results of Broek et al.'^ The figure clearly shows that stanene has a 
Dirac point at K similar to the case in graphene. Thus, it is a zero band gap material. 
Our calculations did not include spin-orbit coupling (SOC) effects. It is known that SOC 
opens up the gap by about 0.1 eV at the K point. In either case, with or without SOC, 
the electron density at the Fermi level is very low for stanene. Consequently, free-standing 
stanene cannot superconduct due to a lack of carriers at the Fermi level that can participate 
in Cooper pairing. 

A first attempt towards making stanene a superconductor is to assume phonon mediated 
superconductivity and mimic the way in which its two dimensional counterpart, graphene 
was made superconducting as reported in Refs. 5,7,16. Hence we doped stanene with an 
electron donor, in this case with Lithium and Calcium, respectively, in order to create an 
electron rich Fermi surface. Electrons at the Fermi surface can then form Cooper pairs and 
make the material superconducting. 

We placed Lithium atoms and Calcium atoms in turn on stanene as shown in Fig. 1 
and optimized the resulting structures but with the hxed lattice constant (a = 4.64 A) 
obtained for free-standing stanene. The optimized structure of Lithium (Calcium) doped 
stanene reveals that the structure is stable when Lithium (Calcium) atoms were placed at 
1.75 (2.24) A above the stanene plane. No imaginary frequencies were seen in the computed 
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phonon dispersion curves (see Fig. 5). 

The electronic band structures of Lithium-doped and Calcium-doped stanene are shown 
in Fig. 4. From the figure, it is clear that carriers are now available near the Fermi surface 
in contrast to the situation with free-standing stanene (Fig. 2). 

A. Eliashberg spectral function and electron-phonon coupling constant 

The Eliashberg spectral function a‘^F{u) and the electron-phonon coupling constant A (a;) 
are plotted as functions of the phonon energy ca in Fig. 6 . Superconductivity in these 
doped stanene systems is clearly in the intermediate coupling regime (i.e., A ~ 1). The 
contributions to the EPC A come from three main ranges of phonon energies: the low energy 
phonons with energies from 0 to ~ 60 cm“^, the intermediate energy phonons (~ 80 — 150 K) 
and the ‘high’ energy regime (^ 150 K). For phonons with zero wavevector (q = 0), it 
is straightforward to determine which phonon modes correspond to each of these regimes. 
However, for those with q 7 ^ 0, it is difficult to unravel the phonon modes that contribute the 
most to a given phonon energy. Thus, it is difficult to determine the modes that enhance A 
the most since all the wavevectors in the hrst Brillouin Zone contribute to the sum in Eq. 4. 
The interaction between the electrons and the low energy phonons contribute more than 
70% to the EPC. The remaining ~ 30% comes mainly from the electronic coupling with the 
intermediate energy phonons. In general, we expect that out-of-plane vibrational modes of 
stanene couples with the tt* and interlayer states to produce superconductivity in the doped 
system. 


B. Transition Temperature 

Table II shows the results obtained from the calculations. We adopted a value of /i* = 0.13 
for the Coulombic repulsion parameter as this value reproduced, in a previous ab initio cal¬ 
culation,^^ the experimental Tc of bulk white Tin. Our results follow the trend of transition 
temperatures observed by Profeta et al.'^ with regard to the ad atoms on a graphene sheet. 
That is, the ad atom that gives the highest is the one whose height above the plane is 
smallest. From Table II, we see that the height 6 of the ad atom above the stanene plane 
is smaller for Lithium doped stanene LiSng than for its Calcium doped counterpart CaSng 
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and, as expected, Tc for the the former (1.3 K) is higher than that of the latter (0.7 K). One 
may argne that the valne of the Conlomb repnlsion parameter nsed in the bnlk case might 
be different from that for the 2D system. In Fig. 7 we show the valnes of Tc determined for 
different valnes of fi* in the range from 0.1 to 0.2 as snggested in Ref. 18. The largest valne 
of Tc obtained was 1.7 K (1.0 K) for LiSng (CaSng). These valnes are signihcantly smaller 
than that predicted for Lithinm-doped graphene (8.1 K), bnt agree with the low valne 
(1.4 K) for Calcinm-doped graphene.'^ 

In conclusion, we reported our results on phonon-mediated superconductivity in Lithium 
and Calcium doped stanene. We presented the electronic band structure, predicted the 
Eliashberg spectral function a‘^F{uj), and showed the contribution of the phonons to the 
total electron phonon coupling coefficient A. Our results reveal that superconductivity can 
be induced in stanene by doping it with alkaline and alkaline earth atoms. We found that 
Tc for Lithium doped stanene should lie in the range from 0.48 to 1.70 K while that of the 
Calcium doped should he somewhere from 0.15 to 1.04 K. These values of Tc are small, 
smaller than in the case for bulk /3-Tin and Lithium-doped graphene. Further work should 
seek for methods to increase Tc either by straining stanene or/and by doping with other ad 
atoms. 
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Table I: Optimized lattice parameter of stanene: a is the lattice constant, d is the out-of-plane 
buckled height, and agn-sn is the Sn — Sn bond length; all distances are in Angstroms 



a 

d 

^sn—sn 

This work 

4.64 

0.92 

2.84 

Broek et al. 

4.62 

0.92 

2.82 

2 

Shou-Cheng et al. 

4.65 

- 

- 
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Table II: Summary of results: 6 (in A) is the height of the ad atom above the stanene plane, ujiog (in 
K) is the logarithmic average of the phonon energy, A is the electron-phonon coupling coefficient, 
and Tc (in K) is the transition temperature. 



5 

^log 

A 

T 

LiSng 

1.75 

60.97 

0.65 

1.3 

CaSng 

2.24 

62.62 

0.54 

0.7 
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Figure 1: Crystal structure of doped stanene. The yellow spheres represent tin atoms while the 
red spheres are the dopants. Both the top (a) and the side views (b) are shown above and below, 
respectively. 
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Figure 2: (a), is the electronic band structure and, (b), is the electronic density of states 

(EDOS) of stanene, Ep is the Fermi energy, which we set as the reference level, and 
F(0, 0,0), M( 2 ,0, 0), and iF(|, 0) are special symmetry k points. 
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Figure 3: Phonon dispersion (left) and the phonon density of states (right) for buckled pristine 
stanene. The phonon dispersion is plotted along P, M if and F high symmetry k points of a 2D 
honeycomb crystal 
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Figure 4: Electronic Band Structures and the Total Density of States. The top panel is for LiSng 
and the bottom one for CaSng. In each case, the electronic band structure is on the left while the 
total electronic density of states (EDOS) is on the right hand side. Ep is the Fermi energy, which 
we set as the reference level with a horizontal line. The band structure is plotted along the k-path 
connecting the T, M Kand F high symmetry k points of a 2D honeycomb crystal 
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Figure 5: Phonon dispersion (left) and the phonon density of states (right) for Lithium-doped 
stanene, LiSng (top) and for Calcium-doped stanene, CaSng (bottom). The phonon dispersion are 
plotted along F, M iLand F high symmetry k points of a 2D honeycomb crystal 
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Figure 6: The Eliashberg spectral function a'^F{io) (solid lines) and the electron-phonon coupling 
constant X{uj) (dotted lines) for Lithium-doped and Calcium-doped stanene. 
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Figure 7: Transition Temperature Tc as a function of 





